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Abstract 

Background: MicroRNAs (miRNAs) play a pivotal role in coordinating messenger RNA (mRNA) transcription and 
stability in almost all known biological processes, including the development of the central nervous system. Despite 
our broad understanding of their involvement, we still have a very sparse understanding of specifically how miRNA 
contribute to the strict regional and temporal regulation of brain development. Accordingly, in the current study 
we have examined the contribution of miRNA in the developing rat telencephalon and mesencephalon from just 
after neural tube closure till birth using a genome-wide microarray strategy. 

Results: We identified temporally distinct expression patterns in both the telencephalon and mesencephalon for 
both miRNAs and their target genes. We demonstrate direct miRNA targeting of several genes involved with the 
migration, differentiation and maturation of neurons. 

Conclusions: Our findings suggest that miRNA have significant implications for the development of neural 
structure and support important mechanisms that if disrupted, may contribute to or drive neurodevelopmental 
disorders. 
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Background 

The developing cerebral cortex is highly compartmental- 
ized and involves precise regulation of cell proliferation, 
migration and differentiation [1,2]. The behaviour of cells 
undergoing these neurodevelopmental processes is gov- 
erned by spatiotemporal changes in gene expression [3], 
Translation and intracellular traffic of these transcripts is 
influenced by small non-coding miRNA that function as 
molecular guide sequences for ribonucleoprotein com- 
plexes involved in targeted mRNA degradation or transla- 
tional repression (reviewed [4-7]). By regulating mRNAs 
post-transcriptional fate, miRNA subsequently play a sig- 
nificant role in organizing complex patterns of gene activ- 
ity. Thus these small RNA molecules also display distinct 
expression profiles [8-10], and exhibit both temporal and 
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spatial specificity [8,11-13]. For example, miR-9 has 
been demonstrated to play a role in controlling neuro- 
genesis timing via targeting of progenitor-promoting 
and cell-cycle exit-promoting genes [14]. Accordingly, 
its expression has been shown to increase during neu- 
rodevelopment, but not in the postnatal period [11], 
and to be regionally restricted to telencephalic, di- 
encephalic and tectal periventricular proliferative zones 
[12]. 

Analysis of miRNA and their target mRNA expression 
in neuronal cell cultures supports a regulatory role for 
miRNAs in neuronal development [15-18]. Consistent 
with this role, a number of neural-specific or neural- 
enriched miRNAs have been identified [19] and are 
thought to be involved in regulating early brain develop- 
ment (reviewed [20]). For example, miR-134 and miR-137 
regulate neuronal maturation through the modulation of 
their target mRNA. The brain-specific miR-134 regulates 
dendritic spine development by repressing the translation 
of an mRNA encoding a protein kinase, Limkl, that con- 
trols spine development [21], while the brain-enriched 
miR-137 has a significant role in neuronal maturation and 



© 2014 Hollins et al.; licensee BioMed Central ftd. This is an Open Access article distributed under the terms of the Creative 
Commons Attribution ficense (http://creativecommons.Org/licenses/by/4.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (http://creativecommons.0rg/publicdomain/zero/1.0/) applies to the data made available in this article, 
unless otherwise stated. 



Hollins et al. BMC Genomics 2014, 15:777 
http://www.biomedcentral.com/1471-2164/15/777 



Page 2 of 13 



dendritic morphogenesis via direct targeting of Mibl, a 
ubiquitin ligase important for neurogenesis and neurode- 
velopment [22]. This supports existing evidence for the 
role of miRNAs in neurogenesis [23], morphogenesis [24], 
neuronal cell specification [25] and oligodendrocyte differ- 
entiation and myelination [26] . 

In the current study we wanted to determine if miRNA 
and their predicted targets are altered in accordance 
with their developmental stage and neurodevelopmen- 
tal regionalization. More specifically we compare devel- 
opmental miRNA and gene expression in the caudal 
regions including the mesencephalon and metencepha- 
lon that develop relatively early, with more rostral areas 
such as the telencephalon and diencephalon which de- 
velop later [27], as we hypothesize that miRNA are the 
subject of related spatiotemporal regulation. By matching 
miRNAs with their inversely-regulated mRNA targets, 
we suggest that miRNAs post-transcriptional regulatory 



influence changes throughout neurodevelopment, with 
a pattern of expression in the developing telencephalon 
occurring after that of the mesencephalon in accord- 
ance with their respective cytoarchitectural and synap- 
tic maturity. 

It is critical for the development and use of rodents as 
animal models of neuropsychiatric disorders that we have 
a good understanding of the nature and timing of regional 
differentiation in the rodent brain. Predicting neurodeve- 
lopmental events and recognizing differing windows of 
vulnerability to environmental insults will help develop 
clinically relevant experimental models. 

Results 

Temporal dynamics of miRNA expression in the 
telencephalon and mesencephalon between El 2 and PO 

miRNA expression was analysed in the telencephalon 
and mesencephalon of embryonic (E) (E12, E14, E15, 
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Figure 1 miRNA expression profiles in the developing rat brain. Intensities of differentially expressed miRNA microarray probes across 
development in the telencephalon (A) and mesencephalon (B). Each colored line represents one probe. Normalized signal intensities range from 
2.3 to 14.3 as indicated by color legend. Color legend is relative to E12 such that probes with low expression at E12 (signal intensity of 2.3) are in 
blue while those with high expression at E12 (signal intensity of 14.3) are in red. (C) Scatter plots illustrating miRNA expression in the 
mesencephalon and telencephalon at each timepoint. Solid line = line of best fit; R 2 = Pearson correlation coefficient. 
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E18, E19) and newborn (postnatal (P) day 0) rats using 
GeneChip miRNA array matrices (Affymetrix). Data ana- 
lysis revealed 303 miRNA in the telencephalon and 266 in 
the mesencephalon with differential expression during 
brain development. miRNA expression displayed greater 
variation in the telencephalon (Figure 1A) compared to 
the mesencephalon (Figure IB) throughout development 
(Figure 1C), with 97 miRNA in the telencephalon having 
substantive temporal variance (SD > 1.3) across develop- 
ment as compared to 38 miRNA in the mesencephalon. 
There were also a number of miRNA with altered expres- 
sion specific to the telencephalon. This occurred at E12 
(32 miRNA) and E19 (26 miRNA) only. In contrast the 
mesencephalon had no region-specific miRNA expression 
at any time during the developmental timepoints exam- 
ined. Furthermore we also observed stage-associated ex- 
pression in the telencephalon-specific miRNA, with 13 
and 7 expressed only at E12 or E19 respectively (Add- 
itional file 1). To gain an appreciation of the biological im- 
plications of these telencephalon-specific miRNA during 
early telencephalon development, putative target genes of 
the 32 miRNA specific to the telencephalon at E12 and 



the 26 miRNA specific to the telencephalon at E19 were 
assessed for Kyoto Encyclopedia of Genes and Genomes 
(KEGG) pathway enrichment using the Gene Annotation 
Tool to Help Explain Relationships (GATHER) online 
database. For the telencephalon-specific miRNA at E12 this 
revealed a number of significantly enriched pathways with 
relevance to neural connectivity and synaptic plasticity, 
such as focal adhesion; adherens junction; oxidative phos- 
phorylation; regulation of actin cytoskeleton and ECM- 
receptor signalling. Similarly, the telencephalon-specific 
miRNA at E19 were predicted to be enriched in the 
Mitogen-activated protein (MAP) kinase signalling path- 
way and Wingless/int (Wnt) signalling pathway (Additional 
file 2). Further analysis using Ingenuity Pathway Analysis 
(IPA) software (Ingenuity Systems) revealed roles in im- 
portant neurodevelopmental processes including: the dif- 
ferentiation of dopamine neurons; neuron and neuroglia 
differentiation; myelination; synaptogenesis; and axonogen- 
esis (Figure 2). 

To determine those miRNA that were differentially 
regulated with respect to developmental age, all miRNA 
were ranked based on their temporal expression variance 
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Figure 2 Computational gene interaction predictions: network of E12 telencephalon-specific miRNA involved in the Neurotransmitters 
and other nervous system signalling pathways. A network of putative gene targets of the 32 telencephalon-specific miRNA at E1 2 was con- 
structed with the Ingenuity Systems Pathway Analysis (IPA) software. Solid grey lines specify direct relationships whereas dotted grey lines in- 
dicate indirect interactions. Genes involved in select neurodevelopmental processes are indicated with solid color lines. 
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through E12 to PO. miRNA undergoing substantial 
changes throughout development (SD > 1.3) were selected 
for downstream analysis (Additional file 3). The 97 
miRNA identified were found to have significandy lower 
expression (F = 5.271, n = 12, p < 0.05) in the 



telencephalon compared to the mesencephalon at embry- 
onic day 14 (Figure 3A). Supervised hierarchical clustering 
of these miRNA was characterized by distinct expression 
clusters, two of which showed significantly lower expres- 
sion (p < 0.001) in the telencephalon compared to the 
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Figure 3 miRNA expression in the developing rat brain. (A) Expression of the 97 miRNA showing strong developmental regulation in the 
telencephalon and mesencephalon throughout development (for display purposes expression was derived by calculating the geomean of signal 
intensity values for all dynamic miRNAs at each developmental timepoint). All data were analysed for significance by one-way ANOVA, followed by 
Tukey-Kramer's post hoc t-testing (* = p < 0.05). (B & C) Relative expression levels for the 97 miRNA showing strong development regulation were 
hierarchically clustered (supervised, correlation uncentred, average linkage; Cluster 3.0). Colors indicate relative signal intensities. (Java Treeview; http:// 
jtreeview.sourceforge.net). Two major clusters are shown in six columns corresponding to the six time points. Comparison of expression of 
telencephalon and mesencephalon miRNA for each cluster is shown adjacent to the heatmap. All data were analysed for significance by one-way 
ANOVA, followed by Tukey-Kramer's post hoc f-test. *** = p < 0.001 . 
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mesencephalon at embryonic day 12 (Figure 3B & C). 
Prominent among this group were miRNAs that are 
enriched in neurons and associated with developmental 
regulation (let-7i, let-7f, let-7b, miR-98), cell cycle regula- 
tion {miR-137, mIR-128) and neural activity (miR-132, 
miR-212). Bioinformatic analysis of these 97 miRNA 
through IPA (Ingenuity Systems) identified 9652 putative 
miRNA target genes. Gene ontology (GO) function and 
KEGG pathway enrichments were performed by mapping 
the predicted target genes from the developmentally regu- 
lated miRNAs using the GATHER online database. For 
these miRNAs and their predicted targets, 22 GO func- 
tions and 3 KEGG pathways were significant (Bayes fac- 
tor >6). The top ontologies included: regulation of cell 
cycle; development; morphogenesis; regulation of tran- 
scription; and neurogenesis (Additional file 4). 

In summary, while we observed substantial overlap be- 
tween miRNA expressed in the telencephalon and mesen- 
cephalon, there was a number of functionally significant 
miRNA expressed exclusively in the telencephalon. From 
longitudinal analysis in each brain region, it was also ap- 
parent that telencephalon miRNA expression was subject 
to more temporal variation compared to the mesenceph- 
alon. Genes predicted to respond to this developmental 
variation in telencephalon cortical miRNA expression 
were enriched in pathways associated with neural develop- 
ment and function. 

Developmental regulation of mRNA expression in the 
telencephalon and mesencephalon between E12 and PO 

mRNA expression was analysed in the telencephalon 
and mesencephalon of embryonic (E12, E14, E15, E18, 
E19) and newborn (PO) rats using RatRef-12 Expression 
BeadChips (Illumina). Data analysis revealed a total of 
16,691 genes with differential expression during neural 
development (E12-P0), comprised of 14,371 genes in the 
telencephalon and 13,906 genes in the mesencephalon. 
To determine those genes that were differentially regu- 
lated with respect to developmental age, all mRNA were 
ranked based on their temporal expression variance be- 
tween E12 to PO. Those that demonstrated the greatest 
variation throughout development (SD >3.2) were selected 
for bioinformatic evaluation (Additional file 5). These 100 



genes with high temporal expression dynamics were fur- 
ther assessed for their involvement in biological pathways. 
GO function and KEGG pathway enrichment of these 
genes revealed 43 significant (p < 0.05) GO functions and 
1 KEGG pathway (Additional file 6). Network analysis 
using IPA (Ingenuity Systems) for predicting potential in- 
teractions between the differentially expressed genes gen- 
erated three networks with the highest IPA analysis scores 
(47, 25 and 21) including Cell-To-Cell Signalling and 
Interaction; Nervous System Development and Function; 
Cellular Growth and Proliferation; and Cellular Develop- 
ment (Table 1). While these biological processes accorded 
well with those predicted by the developmentally regu- 
lated miRNA, we wanted to explore the possibility that 
these 100 dynamically expressed genes were most likely 
responding to changes in miRNA expression and driving 
developmentally significant changes in cell morphology 
and physiology, by bioinformatically associating the puta- 
tive interactions with observed changes in predicted rela- 
tionships or miRNA-mRNA interactions. 

miRNA-mRNA interactions 

Given that miRNAs predominantly function to negatively 
regulate gene expression, we examined the biological sig- 
nificance of inverse expression patterns between differen- 
tially expressed miRNA and differentially expressed genes 
in the telencephalon during development. We used Inge- 
nuity's miRNA Target Filter 8 to analyze all possible 
miRNA and target-gene interactions from the miRNA and 
gene expression array data sets. The results were filtered 
based on the confidence of interaction and an inverse 
miRNA to target mRNA relationship. These gene sets 
were then functionally annotated using the Database for 
Annotation, Visualization, and Integrated Discovery bio- 
informatics resources 6.7 (DAVID) online database, re- 
vealing significant over-representation in neural pathways 
throughout development. This included axon guidance, 
neuroactive ligand-receptor interaction, MAP kinase sig- 
nalling, Wnt signalling, long-term depression and calcium 
signalling (Additional file 7). 

We further investigated the putative roles of mir-98, mir- 
212 and mir-137 in the development of the telencephalon. 
All three miRNA were determined to be dynamically 



Table 1 Network analysis using IPA 

Associated Network Functions Score 

Cell-To-Cell Signalling and Interaction, Nervous System Development and Function, Cell Death 47 

Cellular Growth and Proliferation, Haematological System Development and Function, Haematopoiesis 25 

Cellular Development, Cell Signalling, Organ Development 21 

Cell Morphology, Cellular Function and Maintenance, Cell Cycle 19 

Psychological Disorders, Cell Death, Gene Expression 13 



Network analysis using IPA (Ingenuity Systems) for predicting potential interactions between the differentially expressed genes potentially regulated by 
differentially expressed miRNA. 
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regulated with more than a three-fold increase in 
telencephalon expression between E12 and E19. These 
neuronally important miRNAs have critical roles in neuro- 
development and are associated with neuropathology (see 
Discussion). Concordantly, their predicted gene targets were 
also highly significant in our study. Validation of microarray 
expression for these three miRNA using quantitative real- 
time reverse transcription PCR (Q-PCR) (Figure 4) con- 
firmed significandy lower expression in the telencephalon 
compared to the mesencephalon at E12. To further investi- 
gate their biological relevance, IPA software was used to 
identify putative differentially expressed target genes of 



these dynamically regulated miRNA. Analysis of this gene 
dataset using the DAVID online database revealed signifi- 
cant (p < 0.05) over-representation of neurodevelopmental 
functions in the mammalian forebrain cortex, including 
axonogenesis; neuron projection development; synaptic 
transmission; forebrain development; regionalization; and 
learning and behaviour (Additional file 8). 

Target gene silencing in vitro 

Our expression analysis of miRNA and their cognate 
mRNA targets suggested that there are a tremendous 
number of interactions that underlie the spatiotemporal 
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Figure 4 Select miRNA expression in the developing rat brain. Total RNA was isolated from rat telencephalon and mesencephalon at E12, E14, 
E15, E18, E19 and PO, amplified, labelled and hybridised to Affymetrix Genechip miRNA 1 .0 microarrays. Data were normalised using cubic spline and 
analysed using GeneSpring GX1 1. Figures show the normalised signal intensity values for miR-98 (A) and miR-212 (D) and miR-137 (G) at each 
developmental stage in the mesencephalon and telencephalon (*** = p < 0.001). Findings were validated by Q-PCR and data expressed as a 
fold-change. Mean value of E12 was set to 1-fold induction, and mean values of all other ages were related to E12. MiRNA expression as determined 
by microarray appears as black bars and Q-PCR results appear as white bars. For Q-PCR, relative miRNA expression in the telencephalon was determined 
by the difference between their individual cycle threshold (Ct) value and that produced in the same sample for the U6 snRNA (ACt). Error bars 
represent standard deviation of the means. Microarray expression data was validated by RT-PCR for miR-98 (B), miR-212 (E) and miR-137 (H) with 
correlation coefficients of 0.64 (C), 0.62 (F) and 0.69 (I) respectively. 
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segmentation in mammalian brain development. However, 
to further support the functional association between 
neural development-associated target genes and altered ex- 
pression in this group of miRNAs, the respective miRNA 
recognition elements (MRE) from a selection of six func- 
tionally significant target genes, Clasp2, RablS, Cblll, 
Wnt7a, Gpr88 and Plxnb2, were cloned into the 3' un- 
translated region (UTR) of a luciferase reporter gene 
construct and co-transfected into a recipient cell line 
with synthetic miRNA mimics corresponding to those 
altered in vivo. The degree of reporter gene activity and 
the impact of miRNAs were then ascertained by meas- 
uring the relative luciferase activity (Figure 5). MiRNAs 
interact with their respective MREs in the 3'UTR of tran- 
scripts and each gene may contain binding sites for mul- 
tiple miRNA. Therefore observing miRNA-mediated effects 
using a reporter gene construct allows subtle changes 
in gene expression to be observed. Most of the reporter 
constructs behaved in accordance with expectation and 
were repressed 8-25% (p < 0.05) in the presence of syn- 
thetic miRNA in comparison to a scrambled-duplex co- 
transfected control. These data indicate that transfected 
miRNAs bound to their target MREs and repressed the 
expression of luciferase, with the most responsive tar- 
gets derived from the 3'-UTR of Wnt7a and Gpr88. 
These were both targeted by miR-137, which appeared 
to have the greatest overall effect, whereas miR-212 had 
the least effect on these target gene constructs. The 
only exception was the Plxnb2 construct, which dis- 
played increased expression in the presence of synthetic 
miR-137. Collectively, these reporter assays supported 
the existence of relationships between genes reported 
to be associated with neural development and develop- 
mentally regulated miRNAs. 



Discussion 

miRNA and post-transcriptional regulation of 
neurodevelopment 

The trajectory of human brain development closely mir- 
rors that seen in the rodent [28] and as such the rat of- 
fers an excellent model to investigate the coordination 
and timing of miRNA expression across fetal develop- 
ment. In the current study we examined the expression 
profiles of miRNA within the mesencephalon and telen- 
cephalon from E12 until birth. We hypothesized that the 
developing telencephalon would exhibit a significantly 
different temporal miRNA profile in relation to the mes- 
encephalon and this was supported by the observation 
that 87% of expressed miRNA underwent significant 
change across development, with clear differences in 
patterning between the two regions. The highest level of 
change in miRNA expression occurred early in develop- 
ment (E12-E15) and plateaued thereafter, with differ- 
ences between the telencephalon and mesencephalon 
decreasing noticeably during that time period. In accord- 
ance with our hypothesis, telencephalon miRNA expres- 
sion was significantly lower than the mesencephalon at 
E12 and consistent with the delayed development of this 
region. We also observed 32 miRNA with predicted 
function in neurodevelopmental processes, exclusively 
expressed in the telencephalon during early brain devel- 
opment (El 2). 

We then employed a bioinformatics approach to further 
understand the differential developmental trajectories of 
these two brain regions. Temporally regulated miRNAs 
and genes were examined using a combination of miRNA 
and gene expression profiling and computational target 
prediction. Hierarchical clustering of miRNA during tel- 
encephalon development revealed two distinct clusters of 




□ miR-98 
■ miR-212 

□ miR-137 
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Wnt7a Gpr88 Plxnb2 Rab15 



Figure 5 Regulation of neuronal development-associated reporter gene constructs by miRNA. Chart shows relative activity of reporter 
gene constructs (x-axis) in response to co-transfected miRNAs. Relative luciferase activity for each reporter/miRNA combination was expressed as 
a percentage of the response to scrambled controls. Bars are mean±SEM. (** = p<0.01, * = p<0.05). 
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miRNA. These clusters had significantly higher expression 
in the mesencephalon at E12 whilst expression in the tel- 
encephalon remained significantly lower (Figure 3B & C). 
The miRNA contained in these clusters are predicted to 
target temporally-regulated genes involved in brain signal- 
ling pathways important for neurodevelopment and neur- 
onal functioning including: dopamine receptor signalling; 
GABA receptor signalling; glutamate receptor signalling 
and neuroactive ligand-receptor interactions. As expres- 
sion of the neuroactive ligand-receptor interactions path- 
way is demonstrated to correlate with synaptic maturity 
[29], this data is consistent both with the telencephalon 
being less differentiated at this time point and more im- 
portantly implicates a critical role for miRNA in this 
process. This also indicates that at such an early time- 
point, the developing telencephalon may be largely un- 
affected by environmental insults compared to the 
mesencephalon. This may be particularly relevant for cer- 
tain developmental animal models of disease. For instance 
the widely used maternal immune activation model of 
schizophrenia employs the artificial viral construct poly I: 
C which induces a stricdy time-limited inflammatory re- 
sponse of approximately two days in the maternal animal 
[30] . This may be why when it is commonly administered 
at E9 in mice (equivalent to E10 in a rat), the downstream 
effects may be restricted to the developing mesencephalon 
rather than later maturing regions such as the telenceph- 
alon [31]. 

miRNA-mRNA interactions 

While there are some exceptions, miRNA function is 
usually associated with repression of target gene expres- 
sion, such that when a miRNA is expressed at high 
levels, its target mRNA is reduced [32-34]. Indeed, sev- 
eral reports have identified inverse expression patterns 
between differentially expressed miRNA and differen- 
tially expressed mRNA during development [17,18]. In 
this study we found the putative down-regulated gene 
targets of up-regulated miRNA in the telencephalon be- 
tween E12 and E14 are predicted to be involved in many 
important brain signalling pathways including axon 
guidance, neuroactive ligand-receptor interaction, cal- 
cium signalling and long-term depression. 

Altered expression of miRs-98, -212 and -137 was 
particularly noteworthy, as these three neuronally im- 
portant miRNAs were much more abundant in develop- 
ing mesencephalon than telencephalon at the earliest 
embryonic age showing more than a three-fold increase 
in telencephalon expression between E12 and E19. All 
three miRNA have critical roles in neurodevelopment 
and are associated with neurological disorders such as 
Alzheimer's disease and schizophrenia. We confirmed a 
functional link between these three miRNA and suppres- 
sion of expressed target genes by cloning the MREs from 



Clasp2, RablS, Cblll, Wnt7a, Gpr88 and Plxnb2 down- 
stream of a luciferase reporter gene. These reporter con- 
structs were then co-transfected with synthetic analogues 
of their cognate miRNA into a recipient cell line, such that 
the extent of reporter gene expression and therefore the 
influence of miRNA could be ascertained through the 
measurement of relative luciferase activity (Figure 5). 
Most of these constructs were significandy repressed in 
the presence of synthetic miRNA. In particular, regulation 
of 3'-UTR elements from Wnt7a (miR-137) and Gpr88 
(miR-137) was particularly strong and, along with Clasp2 
{miR-98) and RablS (miR-212) provides support for the 
role of post-transcriptional regulation in axon and den- 
dritic growth, maturation and function [35-38]. Path- 
way analysis of these target genes indicated that these 
miRNA may exert significant influence during neurode- 
velopment through the regulation of pathways including: 
axon guidance; neuroactive ligand-receptor interaction; 
Wnt; Erythroblastic Leukemia Viral Oncogene Homolog 
(ErbB); and MAP kinase signalling. 

In the telencephalon, expression levels of miR-98 in- 
creased over four-fold between E12 and E19, a similar 
finding to that of Sempere and colleagues [15], who 
found mir-98 expression levels progressively accumu- 
lated upon neuronal differentiation in both mouse and 
human embryonic carcinoma cells. Experiments by 
Thomson and colleagues [39] similarly found levels of 
miR-98 were significantly increased in mice at E14.5 
compared to levels in embryonic stem cells. Mir-98 is a 
member of the let-7 family, highly conserved across spe- 
cies in sequence and function and involved in the devel- 
opmental timing of cell fates (reviewed [40]). In support 
of a role for miR-98 in telencephalon development, mir- 
98 overexpression was found to significantly decrease 
the activity of a luciferase reporter gene fused to the 
cytoplasmic linker associated protein 2 {Clasp2) 3'-UTR. 
Clasp2 plays a crucial role in chromosome segregation 
during mitosis [41] and is a key regulator of neuronal 
morphogenesis, neuronal polarity and synapse formation 
and activity [37]. 

Similarly, the results suggest a role for miR-212 in 
brain development with expression levels increased 
more than four-fold in the telencephalon between E12 
and birth. MiR-212 arises from the miR-212/132 cluster, 
which is highly conserved in vertebrates. MiR-212 and 
miR-132 are tandem miRNAs with similar target specifi- 
city due to their identical seed sequences. In the brain, 
miR-212 and miR-132 play vital roles in the formation 
and plasticity of neuronal connections, and long-term 
synapse activation (reviewed [42,43]). Dysregulation of 
miR-212 is associated with drug addiction [44], neuro- 
pathological disorders [43] and schizophrenia [45]. Our 
luciferase reporter assays support miR-212 involvement 
in neurodevelopment via regulation of Rabl5, member 
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RAS oncogene family {RablS) and Cbl proto-oncogene- 
like 1, E3 ubiquitin protein ligase {Cblll). Rabl5 belongs 
to the Rab family of low molecular weight (LMW) G TP- 
binding proteins, crucial in ensuring the spatiotemporal 
regulation of membrane traffic [46]. Rab 15 is enriched 
in neural tissue and plays a role in regulating neuro- 
transmitter release [47]. Cblll is a RING finger type E3 
ubiquitin ligase and an important regulator of cell prolif- 
eration [48]. In Drosophila, Cblll is required for proper 
function of cell surface proteins and for endoderm and 
mesoderm morphogenesis [49]. Consistent with the in- 
crease of miR-212 expression between E12 and E19, 
Cblll function is demonstrated to be most vital in the 
early stages of embryogenesis, with its contribution de- 
clining in later stages [49]. 

Another developmentally regulated miRNA is miR- 
137, with expression levels increasing more than four- 
fold in the telencephalon between E12 and E19. Mir-137 
is enriched in neurons and involved in neuronal matur- 
ation, regulation of dendritic development and pheno- 
typic maturation of new neurons [22]. Recent studies 
demonstrate increasing expression of miR-137 during 
neuronal differentiation, leading to reduced neural stem 
cell proliferation, accelerated neural differentiation and 
decreased dendritic development [22,50]. Dysregulation 
of miR-137 has been associated with intellectual disabil- 
ity and Alzheimer's Disease [51,52] and the MIR137 gene 
has also recently found to be highly associated with 
schizophrenia in a genome wide screen and associated 
with cognitive function in patients with the disorder 
[53,54]. This study confirmed the potency of miR-137 to 
down-regulate G-protein coupled receptor 88 (Gpr88) 
and wingless-type MMTV integration site family, mem- 
ber 7 A (Wnt7a) using a reporter gene strategy. Gpr88 is 
a brain-enriched G protein-coupled receptor with a role 
in modulating the striatal dopaminergic system [55] and 
in regulating medium spiny neuron excitability [38]. 
Wnt7a signalling is critical to regulate dendritic spine 
growth and synaptic strength [35], mediating synapse 
density and numbers, and hippocampal network struc- 
ture [56]. During development Wnt signalling is re- 
quired for cell proliferation and differentiation, cell 
polarity generation and embryonic patterning. 

As various brain regions develop at different times 
there is likely to be differing windows of vulnerability to 
environmental insults. Appropriate spatiotemporal ex- 
pression of miRNA is vital during brain development 
and dysregulation during sensitive stages is likely to have 
a profound impact on the coordinated expression of de- 
velopmentally significant genes associated with the 
pathophysiology of neuropsychiatric disorders. This is 
supported by the alterations of miRNA expression ob- 
served in the neuropathology of Alzheimer's disease 
[51], Huntington's disease [57,58], Parkinson's disease 



[59-61] and schizophrenia [62,63]. Understanding the 
role of miRNA in brain ontogeny could identify critical 
windows of developmental significance that are suscep- 
tible to environmental insults which contribute to the 
development of neurobehavioural or neuropsychiatric 
syndromes. 

Conclusions 

This study represents the first global characterization of 
both miRNA and mRNA expression between the telen- 
cephalon and mesencephalon at various stages during de- 
velopment in rodent brain. Comparing the telencephalon 
and mesencephalon affords a molecular basis for under- 
standing the nature and timing of regional differentiation. 
Although these results may have been hypothesised based 
on current evidence suggesting that the mesencephalon 
develops earlier than the telencephalon, the present data 
now provide direct molecular support that differential re- 
gional and temporal miRNA expression within the brain 
is part of this process. Knowledge of the temporal and 
spatial expression profiles of miRNA and their target 
genes is an important initial step in elucidating their func- 
tions. Together these results show that miRNA expression 
is frequently region- and time-specific with an important 
influence in functionally relevant pathways in the develop- 
ing brain. 

Methods 

RNA preparation and analysis of integrity 

Telencephalon and mesencephalon were dissected from 
brains of embryonic (E) (E12, E14, E15, E18, E19) and 
newborn (PO) rats. Each time point represents a pooled 
litter of rat embryos (n = 8-12). Animals were handled in 
accordance with the University of Queensland's animal 
care and ethics committee approval, in compliance with 
the current Australian Law. Total RNA was isolated from 
the telencephalon and mesencephalon using TRIzol re- 
agent as described previously [64] . RNA concentration, in- 
tegrity and purity were assessed using the Agilent 2100 
Bioanalyser (Agilent Technologies). RNA Integrity Num- 
bers (RIN) were automatically calculated with the pro- 
vided system software [65]. All samples showed RIN 
values superior to 8.5. 

MicroRNA expression profiling 

For miRNA profiling, total RNA (700-1000 ng) was 
labelled using a FlashTag Biotin HSR RNA labelling kit 
according to manufacturer's instructions (Genisphere). 
Labelled RNA was hybridised to the Affymetrix Genechip 
miRNA 1.0 microarrays (miRBase version 11.0) for 
16 hours at 48°C and 60 rpm in a hybridisation oven (Affy- 
metrix) and washed and stained with the Fluidics Station 
450 (Affymetrix). Hybridisation, wash and stain reagents 
were from the Hyb, Wash and Stain kit (Affymetrix). 
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Arrays were scanned with Affymetrix GeneChip Scanner 
3000 7G. Differential miRNA expression was analysed 
using GeneSpring GX 11.0 (Agilent Technologies) and was 
assessed as having a minimum 2-fold difference in either 
direction between two time-points. To determine those 
miRNA that were differentially regulated with respect to 
developmental age, the standard deviation of each miRNA 
across development (E12-P0) was calculated. MiRNA with 
a standard deviation of >1.3 were considered to be devel- 
opmentally regulated. 

mRNA expression profiling 

Total RNA was prepared for microarray using the RNeasy 
Mini Kit (Qiagen) and amplified with the Illumina Total- 
Prep RNA Amplification Kit (Ambion). Microarrays were 
performed using the Illumina Rat-Ref-12_Vl Expression 
BeadChip, following the manufacturer's instructions, and 
scanned using the Illumina BeadArray Reader and BeadS- 
can software (Ambion, Scoresby, VIC, Australia). Bead 
summary data were normalized in GenomeStudio v3 
(Illumina) using Cubic Spline. Differential mRNA ex- 
pression was analysed using GeneSpring GX 11.0 (Agilent 
Technologies) and was assessed as having a minimum 2- 
fold difference in either direction between any two time- 
points. 

Target gene and pathway analyses 

miRNAs that were determined to be differentially 
expressed were further analysed for putative target genes 
under their influence using Ingenuity 8 microRNA target 
filter [Ingenuity 8 Systems (www.ingenuity.com)]. Pathways 
and gene ontology enrichment analysis of the target 
genes was performed using the DAVID (http://david. 
abcc.ncifcrf.gov/) [66,67] and GATHER (http://gath.er. 
genome.duke.edu/) [68] online databases. To further in- 
vestigate biologically relevant targeting, differential 
miRNA and mRNA expression data was analyzed through 
the use of Ingenuity 8 microRNA target filter to generate 
lists of interactions between genes and miRNAs. Results 
were filtered based on a high confidence of interaction 
and an inverse miRNA to target mRNA expression rela- 
tionship. The identified genes were then subjected to 
pathways and gene ontology enrichment analysis. 

Quantitative real-time reverse transcription PCR 

Validation of differentially expressed miRNA was per- 
formed by Q-PCR on a subset of the differentially 
expressed miRNA. Briefly, 500 ng of sample RNA was 
treated with DNase-I (Invitrogen), and multiplex reverse 
transcription was performed with Superscript II reverse 
transcriptase (Invitrogen), a 3 nmol/L mix of miRNA se- 
quence specific primers and primers for U6 small nu- 
clear RNA. The primer sequences for the amplification 
of specific primers were as follows: rno-miR-98 - 5 - 



GTAAAACGACGGCCAGTAACAATA-3 ' (forward), 5'- 
TG AGGTAGTAAGTTG-3 ' (reverse); rno-miR-212 - 5 -T 
AA + CAGTCTCCAG-3 ' (forward), 5 -GTAAAACGACG 
GCCAGTCTACGCG-3 ' (reverse); rno-miR-137 - 5 -TT + 
ATT + GCTTAAGAATA-3' (forward), 5'-GTAAAAC- 
GACGGCCAGTCTACGCG-3 ' (forward); U6 snRNA 
5-TGCTTCGGCAGCACATATAC-3' (forward), 5 -A 
GGGGCCATGCTAATCTTCT-3 ' (reverse). Triplicate re- 
actions were set up in a 96-well format with the epMotion 
5070 automated pipetting system (Eppendorf, Hamburg, 
Germany) and carried out with the Applied Biosystems 
(Foster City, California) 7500 real-time PCR machine. 
Relative quantification was assessed using the formula 2' 
A CT. The delta Ct was calculated by subtracting the Ct of 
the endogenous control (geometric mean of U6 small nu- 
clear RNA (snRNA)) from the Ct of the miRNA. 

Reporter construction 

To generate reporter vectors bearing binding sites for 
the three miRNA examined in detail, miR-98, miR-212 
and miR-137, oligonucleotides encoding target gene 
miRNA recognition elements (MREs) were annealed to 
form Spel and Hindlll restricted overhangs of a ligatable 
cassette compatible with Spel and Hindlll digested 
pMIR-REPORT vector (Ambion, Austin, TX) as de- 
scribed previously [69,70]. Recombinant reporter con- 
structs were amplified in chemicompetent DH5alpha 
Escherichia coli bacteria cells (Bioline, Sydney) and puri- 
fied using endotoxin free minipreps (Promega, Madison, 
WI) before being verified by DNA sequencing (AGRF, 
St. Lucia, QLD, Australia) and quantitated using a Nano- 
Drop 2000 (Thermo Scientific, Delaware, ME, USA). 

Cell culture 

Human embryonic kidney (HEK)-293 cell cultures were 
maintained as confluent monolayers at 37°C with 5% 
CO2 and 90% humidity in Dulbecco's modified Eagles 
medium (DMEM) with 10% (vol/voP 1 ) fetal calf serum, 
20 mM HEPES (4-(2-hydroxyethyl)-l-piperazineethane- 
sulfonic acid), 0.15% (wt/vor 1 ) sodium bicarbonate and 
2 mM L-glutamine. 

Transfection 

For the luciferase assay, HEK-293 cells were seeded into a 
96-well plate at a density of 20,000 cells per well and 
transfected 24 h later using Lipofectamine 2000 (Invitro- 
gen). In each case transfections were performed according 
to the manufacturer's instructions, with lOOnM synthetic 
miRNA (Invitrogen). Validation of predicted target 
genes was accomplished by co-transfecting HEK293 
cells with synthetic miRNA and recombinant firefly lu- 
ciferase reporter gene constructs containing 3'-UTR se- 
quences substituted from the target gene. Each reaction 
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identified using IPA. Only those pathways and ontologies with a Bayes 
factor >6 were considered significant. 

Additional file 3: A table listing the 97 dynamically regulated 
miRNA. Also listed are the standard deviation values for each miRNA in 
both the forebrain and midbrain. 

Additional file 4: A table listing enriched ontologies of putative 
target genes of the 97 dynamically regulated miRNA identified 
using the GATHER online database. Putative gene targets were 
identified using IPA. Only those pathways and ontologies with a Bayes 
factor >6 were considered significant. 

Additional file 5: A table listing the 100 dynamically regulated 
mRNA. Also listed are the standard deviation values for each mRNA in 
both the forebrain and midbrain. 

Additional file 6: A table listing enriched ontologies of the 100 
dynamically regulated miRNA identified using the GATHER online 
database. Only those pathways and ontologies with a Bayes factor >6 
were considered significant. 

Additional file 7: A table listing enriched pathways of differentially 
expressed miRNA targeting differentially expressed RNA identified 
using the DAVID online database. Only those pathways with a p-value 
<0.05 were considered significant. 

Additional file 8: A table listing predicted ontologies enriched with 
target genes of miRs-98, -212 andX -737. Enriched ontologies of 
putative target genes of miR-98, miR-212 and miR-137 were identified 
using the DAVID online database. Putative gene targets were identified 
using an intersection of the miRanda, IPA and TargetScan output, n, 
number of input genes in pathway.%, percentage of total target genes 
analysed that are represented in this pathway. 



was performed in quadruplicate and replicated at least 
once. To validate the putative MREs, the HEK293 cell 
line, which is of embryonic origin and easily transfected, 
was an ideal candidate for the investigation of the func- 
tional regulation of miRNA. Candidate genes were se- 
lected based on: their inverse expression respective to 
their targeting miRNA; their role in neurodevelopment 
and association with neuropathology; a high MRE-binding 
specificity score using the miRanda online software [71] 
{Clasp2 - 142; RablS - 158; Cblll - 161; Wnt7a - 149; 
Gpr88 - 150; Plxnb2 - 161). 

Reporter gene assay 

Twenty-four hours after transfection, cells were har- 
vested in 100 ml Passive Lysis Buffer (Promega, Madi- 
son, WI) and 10 ml lysate was used to measure the 
relative luciferase activity using the dual-luciferase re- 
porter assay kit (Promega, Madison, WI). Reporter gene 
silencing in response to miRNA co-transfection was 
monitored with respect to a pRL-TK control plasmid ex- 
pressing renilla luciferase using a Synergy plate reader 
with dual reagent injector (Biotek). To control for the 
nonspecific effects associated with transfection, the con- 
trols were co-transfected with mutant miRNAs or mu- 
tant anti-miRs (Invitrogen). 

Statistical analysis 

Statistical significance analysis was assessed using one-way 
or two-way analysis of variance (ANOVA) with Benjamini 
Hochberg multiple comparisons testing (Prism software 
version 5.0; GraphPad Software Inc., San Diego, CA). A 
value of p < 0.05 was considered significant. The Tukey- 
HSD test with a level of significance set at p < 0.05 was 
used for post hoc comparisons. Unsupervised hierarchical 
clustering was performed in Cluster (Stanford University, 
Palo Alto, CA, USA) [72] . Data were log transformed and 
median centred by genes. Genes and arrays were clus- 
tered, correlation uncentred, by average linkage clustering 
and visualized through Java Treeview V.1.1.6r2 (http:// 
jtreeview.sourceforge.net) [73]. Statistical comparison of 
luciferase data sets was performed by two-tailed paired 
Student's i-test. 

Additional files 
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Additional file 1: A table listing miRNA expressed in the 
telencephalon but not the mesencephalon at El 2 and at El 9. Also 
listed are miRNA that are specific to the telencephalon at E12 only and at 
E19 only. 

Additional file 2: A table listing putative pathways regulated by 

the El 2 and El 9 telencephalon-specific miRNA. Enriched pathways of 
putative target genes of the 32 miRNA specific to the telencephalon at 
E12 and the 26 miRNA specific to the telencephalon at E19 were 
identified using the GATHER online database. Putative gene targets were 
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